Thermal Casimir effect for Drude metals in the plane-sphere geometry 
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We compute the Casimir interaction between a plane and a sphere, the configuration employed 
in the most precise experiments. The scattering formula is developed by taking a suitably chosen 
plane- wave and multipole basis and is valid for arbitrary values of the sphere radius, inter-plate dis- 
tance, temperature and arbitrary dielectric functions for both sphere and plate. Our analytical and 
numerical results for metallic surfaces show a non-trivial interplay between the effects of curvature, 
temperature, finite conductivity and dissipation. 
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I. INTRODUCTION 

Measuring the Casimir force [l|, |2| has been the aim 
of an increasing number of experiments in the last years 
|3l4l7j| . The comparison of these measurements with pre- 
dictions from quantum electrodynamics theory have been 
applied to put constraints on hypothetical new forces pre- 
dicted by unification models [lg, Gil • Accurate theoret- 
ical computations, accounting for a modeling of experi- 
mental conditions, are sorely needed for all comparisons 
to be reliable [2(1 [HI- 

The Casimir force is indeed very sensitive to experi- 
mental conditions. The effect of finite conductivity [22[ 
plays an essential role in the accurate determination of 
the force while the contribution of thermal fluctuations 
give s rise to a remarkable interplay with the latter effect 
[23l - l28l | - In the calculations performed for the geome- 
try of two parallel plates, the Casimir force computed 
within the dissipative Drude model turns out to be a 
factor of 2 smaller than the result obtained within the 
lossless plasma model. As a consequence, the plasma 
theory of the Casimir effect cannot be obtained from the 
Drude one by taking the corresponding relaxation con- 
stant to zero. The current experimental results [Tl| do 
not explore the region where the calculations for lossy 
and lossless models give large differences but their pre- 
cision is sufficient to favor the plasma over the Drude 
theory, in spite of the dissipative nature of the metallic 
plates used in the experiments. 

Now, the most precise experiments are performed with 
a spherical metallic surface in front of a plane surface. 
The force is usually derived from the force evaluated 
in the parallel-plate geometry with subsequent averag- 
ing over the local separation distances. This Proximity 
Force Approximation (PFA) '29>] is expected to provide 
an accurate description in the limit of large sphere radii 
(see [30j, |3l| for derivations with perfect mirrors at zero 
temperature). Even if spheres used in the experiments 
are much larger than the typical distance between them, 
it remains necessary to master the beyond-PFA geometry 
correction [32j even for large spheres in order to match 



the experimental accuracy level. 

There is no reason why the thermal, finite-conductivity 
and beyond-PFA corrections could be expected to be in- 
dependent. Therefore, an accurate description of the ex- 
perimental conditions has to take these effects into ac- 
count simultaneously within a single theoretical model. 
In this paper, we develop the beyond-PFA scattering ap- 
proach [20( for the plane-sphere geometry at finite tem- 
perature, with material properties described by either the 
perfect reflector, plasma or Drude models. We show that 
the interplay between temperature and material proper- 
ties is drastically affected when the parallel-plate config- 
uration is replaced by the plane-sphere geometry. The 
results obtained from the Drude and plasma models are 
generally closer to each other than in the parallel-plate 
geometry. In particular, the factor of 2 between the two 
models reached at the limit of long distances between 
parallel plates is reduced to 3/2 for a plate and a sphere, 
and even less if small spheres are considered. Finally, 
we find that PFA underestimates the thermal contribu- 
tion to the Casimir force for the Drude model at short 
distances, whereas it overestimates it at all distances for 
the perfect reflector and plasma models. 

We start from the scattering theory of Casimir inter- 
action [20], which allows one to consider non-trivial ge- 
ometries at finite temperatures together with a realistic 
description of the material properties. The Casimir free 
energy is written in terms of general reflection operators, 
which describe non-specular diffraction by the material 
surfaces. The resulting formula provides a compact way 
of taking into account the multiple scatterings between 
the interacting bodies [33| • In the particular case of par- 
allel plane surfaces, the reflection operators are diagonal 
in the plane-wave basis, so that the Casimir free energy 
is given in terms of specular reflection coefficients [34J . 
In the plane-sphere geometry, scattering on the sphere 
is easily made explicit in the multipole spherical-waves 
basis [35[ (see also [36[ for a treatment of non-retarded 
interactions). It is also essential to combine the multi- 
pole basis with the plane wave basis |37j ] well adapted 
to the treatment of non-ideal reflection by metals. By 
judiciously employing the two bases at appropriate steps 



of the derivation, the scattering formalism thus allows to 
analyse the inter play between geometry and finite con- 
ductivity at zero [38[ and non-zero temperature [39j j . 

Let us note that several papers have been devoted in 
the recent years to the study of the Casimir effect in non- 
planar geometries. Besides the already cited papers, let 
us refer to the following applications of the scattering 
approach |40l - l47| or of alternative methods 48 55]. The 
first application of the scattering approach to non-trivial 
geometries and non-ideal reflectors can be found in [5fl- 
l58j ]. An historial overview of the various forms of the 
multiple scattering method can be found in [59(. In this 
paper, we analyse in a detailed manner the full interplay 
between finite conductivity, finite temperature and plane- 
sphere geometry, which was only briefly discussed in [39j . 

This paper is organized as follows: our main goals and 
basic notations are presented in Sec. II. In Sec. Ill, 
we develop the scattering approach for the plane-sphere 
geometry and derive the formal results used in Sees. IV 
(perfect reflectors) and V (plasma and Drude models). 
Concluding remarks are presented in Sec. VI. 



II. DEFINITIONS AND GENERAL OUTLOOK 

We consider the Casimir interaction between a metallic 
sphere of radius R and a metallic plate at a distance 
of closest approach L at an arbitrary temperature T, as 
shown in Fig. 1. The center-to-plate distance is 

£ = L + R. 

Experimental results for the Casimir force in the plane- 
sphere geometry are usually compared with PFA-based 
theoretical models [29[ . The spherical surface is assumed 
to be nearly flat over a scale of the order of L, and then 
the total Casimir energy is obtained by adding the con- 
tributions corresponding to different local inter-plate dis- 
tances over the sphere surface. Although the Casimir 
energy is not additive, PEA is usually expected to pro- 
vide an accurate description when R^> L. The resulting 
Casimir force F PFA is 



F 



PFA 



2irR 



J-pp 

~A~' 



(1) 



J-pp is the Casimir free energy for two parallel plates 
of area A. Thus, PFA neglects diffraction due to the 
sphere curvature and provides a direct connection to the 
parallel-plates geometry, the force being proportional to 
the Casimir free energy per unit area calculated for this 
much simpler geometry. 

In this paper, we develop an exact theoretical model 
taking diffraction fully into account. We first compute 
the Casimir free energy J-(L, T) for the plane-sphere ge- 
ometry from the scattering formula [2fJ (Sec. Ill) and 
then derive the force F and entropy S from 
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FIG. 1: Sphere of radius R and a flat plate at a distance L, 
with center-to-plate distance C = L + R. 



The latter is at the core of the ongoing debate about a 
possible violation of the third law of thermodynamics in 
the dissipative Drude model [25|, |27J]. Negative Casimir 
entropy values found for parallel plates have been ex- 
plained in terms of the coupling to a heat bath associated 
to dissipation [28]. But here we find negative entropies 
even for the perfect reflector model, showing that ge- 
ometry itself plays a non-trivial thermodynamical role, 
provided that beyond-PFA diffraction effects are taken 
into account. 

To quantify the deviation from PFA, we define the 
quantity 



F 



PF 



The ratio 



^PFA' 



F(L,T) 
F(L,0) 



(3) 



(4) 



represents the temperature correction at a given separa- 
tion distance and for a given model. In the next sections, 
we calculate pp and -d for the perfect reflector, plasma 
and Drude models for metallic surfaces. 
The Drude dielectric function 



<i0 



m+7) 



(5) 



(2) 



at imaginary frequencies uj — i£ contains two frequency 
scales: the plasma frequency ujp and the relaxation fre- 
quency 7. The plasma dielectric function is obtained from 
([5]) in the lossless limit 7 — > 0. Note, however, that there 
is no continuity in the Casimir force from the Drude to 
the plasma model in this limit [28|. By taking the fur- 
ther limit ujp — > 00 from the plasma model, we recover 
the perfect reflector limit (to be discussed in Sec. IV), 
which corresponds to an infinite dielectric function at all 
frequencies. 

In Sec. V, we provide a detailed comparison between 
plasma and Drude models for the Casimir effect. Since 
the dc conductivity <jq — utp/j diverges in the limit 
7^0, the Drude model is expected, in principle, to 
provide a more realistic description of normal (i.e., non- 
superconducting) metals. However, experimental data 
are surprisingly in better agreement with the plasma 



model when PFA is employed to analyze the plane-sphere 
geometry (llj. Here we show that the results from the 
two models are actually closer than predicted by PFA- 
based theories, which might help to solve this paradox. 
We consider a rich parameter space containing five dif- 
ferent length scales: besides the geometrical scales L and 
R, our problem contains the thermal wavelength (ks is 
the Boltzmann constant) 



At = 



he 



the plasma wavelength Xp — 2t:c/ljp and the wavelength 
corresponding to the relaxation frequency A 7 = 2irc/"f. 
Different orderings of these length scales can in principle 
be considered, leading to various regimes associated to 
non-trivial interplays between temperature, geometry, fi- 
nite conductivity and dissipation. Most experiments are 
performed with gold at room temperature, hence the nu- 
merical results presented in the following sections corre- 
spond to Xt = 7.6 /mi, Xp — 136 nm (plasma and Drude) 
and A 7 /A P = 250 (Drude). 



III. SCATTERING APPROACH IN THE 
PLANE-SPHERE GEOMETRY 

In this section, we apply the scattering approach to the 
plane-sphere geometry (see Fig. 1) at finite temperature. 
The Casimir free energy is written as a sum over the 
Matsubara frequencies £„ (n > 0) [201 ] : 



lixnliBT 



T = k B TY^ log det (1 - M(£ n )) , £„ = 

n 

M(0 = Ks(0e- K M c K P (£)e- K ^ c , 



(6) 



where the primed sum means that the n = term is 
counted for a half. The reflection operators of the sphere, 
7£s(£)i an d the plate, 7£p(£), are evaluated with refer- 
ence points at the sphere center and at its projection on 
the plane, respectively. The operator e~ !C ^ c accounts 
for one-way propagation along the z axis between these 
points, separated by the length C. Thus, the operator 
M(0 represents one round-trip propagation inside the 
open cavity formed by the two surfaces. 

The plane- wave basis |k, ±,p)^ (k = wavevector com- 
ponent parallel to the xy plane, p — TE, TM for po- 
larization and +/— for upwards/downwards propagation 
direction) is well adapted to the description of the prop- 
agation operator e~ lc ^ c , which is diagonal in this basis 
with matrix elements e~ KC , k — \/£, 2 /c 2 + k 2 represent- 
ing the wave- vector z-component associated to the imag- 
inary frequency £. Reflection on the plane also preserves 
all plane wave quantum numbers but the propagation di- 
rection, and the non zero elements of lZp(£) are given 
by the standard Fresnel specular reflection amplitudes 
r p (k, £) for an homogenous medium. 



On the other hand, the multipole basis \£mP)^, with 
£(£ + 1) and m denoting the usual angular momen- 
tum eigenvalues (with £ = 1,2,..., m = —£,...,£) and 
P = E,M representing electric and magnetic multipoles, 
is well adapted to the spherical symmetry of IZs (£) • By 
rotational symmetry around the z-axis, A4(£) commutes 
with the angular momentum operator J z . Hence .M(£) is 
block diagonal, and each block M ^ m ' (£) (corresponding 
to a given subspace to) yields an independent contribu- 
tion to the Casimir energy. We find its matrix elements 
in the multipole basis after introducing the spectral res- 
olution of the identity operator in the plane wave basis: 



X (m) (0l.2 = 



d 2 k 



J2 (£imP 1 \TZ s (0\K+,p) 

p=TE,TM 

xr p (k,e)e- 2KX (k,-,p|^mP 2 ). 



(27T) 



(7) 



This expression has a simple interpretation when read 
from right to left: a multipole wave \£imP-2) is first de- 
composed into plane waves (coefficients (k, — ^p^mP'i)) 
which propagate towards the plane (factor e~ KC ). Af- 
ter reflection by the plane (specular amplitude r p (k, £)), 
the plane wave components propagate back to the sphere 
(second factor e~ KC ) and are finally scattered into a new 
multipole wave \£\mPi). 

The matrix elements of 1Zs in Q represent the multi- 
pole components of the field scattered by the sphere for a 
given incident plane wave. In Mie scattering calculations, 
one usually assumes that the incident plane wave prop- 
agates along the z-direction [61(. Here this is no longer 
possible, since we have to consider field modes propagat- 
ing simultaneously along all possible directions. It is then 
useful to re-formulate [62j the Mie scattering expressions 
in terms of the matrix elements of finite rotation 63] 



% ,{d) = (£m\e- l9J y\£m'} 



with ml = ±1 accounting for the photon spin. 

The resulting expressions for (£imPi|7vls(£)|k, +,p) 
(see appendix), are proportional to the Mie coefficients 
ai(i£) and bg(i£) [611 ] . which represent the scattering am- 
plitudes for electric and magnetic multipole waves. At 
the imaginary frequency axis, they are written in terms 
of the modified Bessel functions [64| evaluated at the 'size 
parameter' 



£ = 



£R 



as follows: 




ae(i£) = 


2 0) 0) 


2 n 2 s[ — s\ 


bM) = 


(a) (6) 

r, („\ l r l\ 



(8) 
(9) 



» _ 



„(*>) 



» 



M) _ 



4+l/ 2 (<) (1^+1/2(1) -^-i/ 2 (D) 
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with n = \/e representing the sphere refractive index. 

We also derive in the appendix the change-of-basis ma- 
trix elements (k, — ,p\l 2 mP%). They yield, when replaced 
into (|7|), explicit expressions for the matrix elements of 
the round-trip operator M.( m '(£,), which we organize as 
a block matrix: 



M^(F\=[ M(m) ( E > E ) MW(E,M) 

u " l A^ m )(M,£) M {m) {M,M) 



(10) 



Each block is the sum of TE and TM contributions: 
M^(Pi,P 2 ) = T, p M { ™ ) (Pi,P 2 ). For the diagonal 
blocks, we find 
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The matrices A^ m ^ and £?( m ) do neither depend on the 
radius R nor on the refractive index of the sphere (spher- 
ical harmonics Yi m (9, if) [63j calculated at the azimuthal 
angle if — 0): 
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' r ~dfc-[e i ^ (+) )-^.- 1 (^ (+) )](16) 
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±-, KSv /e 2 /c 2 + fc 2 .(17) 



Similar expressions are found for the non-diagonal 



blocks: 

M^(E,M) ei , £2 = iy/%0SCJ£l TE a t M) 
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(7( m ) anc j £)("») are a i so written in terms of spherical 
harmonics and rotation matrices: 
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By inspection of ([T5j) and (fTB]) . it is easy to show 
that the diagonal blocks are invariant when we replace 
to — > —to,, whereas the non-diagonal blocks given by (|18[) 
and p9[) change sign. Thus, the matrix A^ ** is block 
diagonal and det(l — M^(^ n )) does not depend on the 
sign of to, allowing us to write the Casimir free energy 
as a double primed sum (with only non-negative values 
of to): 



T = 2£; B T^^logdct (l - X (m) (£„)) 



(20) 



The zero-temperature limit can also be obtained from 
([20jl by changing the sum over Matsubara frequencies by 
an integral over £ : 



^rV^sl $. 






2tt 



(r->o). (21) 



Numerical considerations. When evaluating the free 
energy from (|2U)l . one needs to truncate the vector space 
at some maximum value ^ max of angular momentum. 
From the localization principle [651 ] , the value of £ max re- 
quired for a given accuracy level is expected to scale with 
the size parameter £ which captures the dependence of 
_/Vf( m ) on the sphere radius according to ©-([9]). Charac- 
teristic frequencies, giving the main contribution to the 
Casimir energy, scale as £ ~ c/L. As a consequence, the 
required ^ max should scale as R/L for intermediate and 
short separation distances. This is verified by the nu- 
merical evaluations discussed below. They are limited to 
! max = 85 at the moment (^ max = 45 for Drude model), 
allowing us to calculate for R/L < 20 (R/L < 10 for 
Drude model), not far from for the experimental range 
[22] R/L > 10 2 which would require a larger £ max . 

Apart from this restriction, we are able to calculate the 
exact Casimir free energy and force for parameters closely 




FIG. 2: Casimir force as a function of distance L at T = 
300K for a sphere radius R = 10 /irn, computed for perfect 
reflectors (dashed) or gold surfaces described by the Drude 
model (solid). The dotted line represents the PFA result for 
perfect reflectors. 



mimicking the experimental conditions. As a first exam- 
ple, we plot in Fig. [2]thc Casimir force as a function of L 
for a sphere of radius R = 10 /im at room temperature. 
The dashed and solid lines correspond to the perfect re- 
flector model (to be discussed in detail in Sec. IV) and 
to the Drude model (Sec. V), respectively. We also show 
the PFA result for perfect reflectors (dotted line). The 
interplay between the effects of temperature, geometry, 
finite conductivity and dissipation is better understood 
when considering the ratios defined in Sec. II, which are 
analyzed in detail in the next sections. At the limit of 
large separation distances, simple analytical results can 
be derived, as discussed below. 

Large- distance limit. When R <C £, the characteristic 
size parameters scale as £ ~ R/C -C 1 and the Mie co- 
efficients at(i£) ~ £, 2e+1 , with the magnetic coefficients 
bt(i£) of the same order or smaller than a^(i£) depending 
on the material properties of the sphere. The resulting 
matrix elements My n '{P, P')e ly e 2 are very small and the 
dominant contribution comes from I = 1 (m = 0,1), 
which corresponds to the dipole contributions: 



logdct(l-X^)«- Yi, M (m) (P,P)i 



(22) 



P=E,M 



The explicit expressions for the electric dipole matrix el- 
ements are obtained from ([TT ]) -([T2" ]) and ([15 ]) - ([IS ]) : 



M {0 \E,E) 1A 
M^(E,E) 1A 



3ai(i£) r°° 
2(C/c) 3 JO K 



3ai(iQ 
4(5/c)3 



Jo d* -rTue- 2KC (23) 
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The magnetic dipole elements are obtained from (j2"3")l - 
flU} by replacing ai(i£) — > bi(i£) and interchanging 
ttm <-> fTE- When the refractive index is finite, its contri- 
bution is negligible and then the dominant contribution 



is the electric dipole one (Rayleigh scattering). To pro- 
ceed further, we need to specify the material properties 
for the sphere and plane and the corresponding Mie and 
Fresnel coefficients appearing in (f2"3" |) - (|2"4")) . In the next 
two sections, we consider perfect reflectors, plasma and 
Drude metals. 



IV. PERFECT REFLECTORS 

Within the perfect reflector model for metals, the di- 
electric permittivity is taken to be infinite at all frequen- 
cies. This simple albeit unphysical model for metals pro- 
vides an accurate description at large separations in the 
ideal zero-temperature case. At finite temperatures, how- 
ever, it is still unclear if the perfect reflector model repro- 
duces the correct long distance regime for real metals. In 
fact, for parallel plates, it predicts a force twice as large 
as the value obtained within the dissipative Drude model 
for metals in the long-distance limit, while we would in 
principle expect the two results to agree in this limit. 

It is thus extremely important to compare the re- 
sults obtained from the different models for the plane- 
sphere geometry. In this section, we start with the per- 
fect reflector model, with Fresnel reflection coefficients 
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-1. The Mie coefficients are obtained 



by taking n > 1 and n£ > 1 in ©-©: 



pcrf 



m = in 



e+ 



! iit+iMi)-tit-iMi) 



^+1/2(0+^-1/2(0 



k( 2 5) 



Z ^+1/2(0 



(26) 



The Fresnel and Mie coefficients written above can also 
be obtained from the plasma model expressions (dis- 
cussed in the next section) by taking the limit Ap — > 0. 

Large- distance limit for perfect reflectors. For £ -C 1 
we may take the power expansion of (|25p - (|2"6"]) : 



( '0 = ~l? + O(f) 5 , bT l( (i0 = I f 3 + O(0 5 . 



, pcrf /,-A _ 



Note that b\ , representing the magnetic dipole contribu- 
tion, is of the same order of the electric dipole coefficient 
a\ for £ <C 1. This property holds whenever n£ ^> 1, 
which is also the case for the plasma model at low fre- 
quencies when Ap -C R. On the other hand, for any finite 
dielectric constant, the magnetic dipole is always of the 
order of the electric quadrupole and thus much smaller 
than the electric dipole contribution in the low- frequency 
limit. We find this so-called Rayleigh scattering regime 
when discussing the long-distance limit within the Drude 
model in the next section. 

We insert ((27)) and the values for the Fresnel coeffi- 
cients into (|2T)|) - (|2"41 and write the resulting expression for 



the free energy in terms of the thermal wavelength At : 
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The evaluation of the sum over Matsubara frequencies in 
is straightforward: 
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_ ^ sinh v + cosh i/(^ 2 + sinh v) 
2 sinh 3 v 



The low temperature approximation is derived by ex- 
panding (f>[y) in powers of v : 
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, A T > C > i?. 



. . (31) 
Note that the zero-temperature limit contained in the 

above expression can also be obtained by replacing the 

sum over n in (1281) by an integral over £ as in (|21[) . 

The first finite-temperature correction in pip does not 
depend on C, so that the temperature correction to the 
Casimir force comes from the next-to-leading-order term, 
proportional to (fc^T) 6 . It is a repulsive contribution 
from thermal photons that makes the net force slightly 
less attractive. We discuss this point further in the con- 
text of the numerical evaluations presented below. 

We can compute the high-temperature limit by taking 
v — > oo in (130[) or by considering the n = Matsubara 
frequency contribution in 
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(32) 



It is also interesting to compute the Casimir entropy 
from (1301): 



5 pe 



3k B R 3 
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(tf>(i/)+i/0'(i/)), £>#. (33) 



Remarkably, this expression yields S p < for v < 1.5, 
that is L < 1.8/xm at T = 300K. In the parallel plates 
geometry, negative entropies have been found for the dis- 
sipative Drude model only [60j ■ Here we find negative en- 
tropies also for perfect reflectors over a wide separation 
distance range, as discussed below. 

Numerical results. In Fig. [3p,, we plot the ratio 
i9 porf [see 0], quantifying the thermal correction to the 
Casimir force, for different sphere radii as a function of 
the separation distance L. We also show the results ob- 
tained by using the PFA formula ([1]) (dashed line) and by 
deriving the analytical long-distance expression (j3"0"|) with 
respect to L after replacing C « L (dotted line). The 
agreement with the analytical formula is as expected bet- 
ter if we keep the variable C in (f30)) . as shown in Fig. [3Jd, 
where we plot $P crf as a function of C. 
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FIG. 3: Thermal Casimir force at T = 300K computed be- 
tween perfectly reflecting sphere and plane, divided by the 
zero temperature force, as function of surface distance L (a) 
and distance-to-center C (b). Solid lines from bottom to top 
correspond to increasing values of sphere radii. The upper 
dashed curve represents the PFA result while the lower dotted 
curve corresponds to the analytical large-distance expression. 
[Colors online]. 



At very short distances L < At we recover, as ex- 
pected, the zero temperature result (d — 1). As the 
distance increases, we find, in most cases, that § de- 
creases below one, reaches a radius-dependent minimum 
and then increases again at long distances. As long as 
R is not too large, the thermal photons provide a repul- 
sive contribution (thus decreasing the magnitude of the 
overall attractive force) over a distance range that be- 
comes wider as R decreases. This range corresponds to 
L < At/7t when R <§C L, as obtained from (|30|) . 

The reduction of the Casimir force is consistent with 
the negative Casimir entropies found from (1331) in the 
limit R <§; L. We plot S as a function of L (Fig. fjh.) 
or C (Fig. Hfc) for different sphere radii. Fig. HJd shows 
that (|33)) provides an accurate description for C/R > 4. 
Negative entropies are found for R as large as 1 /im. 

An additional relevant property can be inferred from 
Fig. [31 the PFA always overestimates the thermal cor- 
rection for perfect reflectors, the overestimation being 
smaller for larger radii (at a given separation distance 
L) as expected. 
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FIG. 4: Casimir entropy at T — 300K computed for perfectly 
reflecting sphere and plane, as a function of surface separation 
distance L (a) and center-to-plate distance C(b). The dashed 
curve corresponds to the analytical asymptotic expression for 
£ 3> R. [Colors online]. 



In Fig.O we plot the beyond-PFA correction p at room 
temperature as a function of L/R, for different values of 
R. At zero temperature, the different curves shown in 
this figure would collapse into a single one |3a, [37| ■ For 
a given ratio L/R, the thermal reduction effect already 
apparent in Fig. [3] is larger for larger radii. The fact 
that pF depends on R for a given L/R is again a clear 
signature of the interplay between thermal and geometry 
effects, that can damage the precision of PFA. In the 
next section, we discuss how this interplay is modified 
when finite conductivity and dissipation are included in 
the model. 



V. PLASMA AND DRUDE METALS 

The plasma model provides the simplest way to take 
the finite conductivity of metals into account. The 
Drude model is a more accurate description of non- 
superconducting metals since it also includes the relax- 
ation of conduction electrons and the associated finite dc 
conductivity. In this section, we take the dielectric con- 
stant e given by (JSJ) (with 7 = in the plasma model) and 
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FIG. 5: Beyond-PFA correction factor computed for perfectly 
reflecting sphere and plane at ambient temperature (solid 
lines) and for asymptotic cases (dashed and dotted lines), as 
a function of surface separation L. [Colors online]. 



derive analytical and numerical results for the Casimir 
free energy and force. 

Large distance high-temperature limit. When C 3> R, 
we take the low-frequency expansion of the Mie coeffi- 
cients and find, for the plasma model [66} , 



plas 



2£ 3 
3 ' 




f 1 l 
2irR 


cotha 
a 



Xi- 



e 3 (34) 



(35) 



Here we also assume that C ^> At (high-temperature 
limit), so that we take only the first Matsubara frequency 
£o = when computing the Casimir free energy from 
(|20l) In the low frequency limit, the Fresnel coefficients 
are given by 7"te ~ -^tm ~ — 1, and then we find, from 
plugging ((Ml) into Eqs. flt^-flM]) : 



^pla 



3hcR 3 
8X T C 3 
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cothc 



a 



C-»\ T ,R. 



(36) 

This result reproduces, as a particular case, the perfectly- 
reflecting limit given by (|3"2"j) when Ap <C R. 

For the Drude model, the TE Fresnel reflection coeffi- 
cient has the well-known low-frequency limit tte ~> 0, 
whereas the TM coefficient behaves as in the plasma 
model: ttm ~ 1- The low- frequency expansion of the 
Mie coefficients are also quite different from the plasma 
case. For the electric dipole coefficient, we find 



,Drud 



2^ 
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a R' 



(37) 



As discussed in Sec. IV, the magnetic dipole b\ is always 
much smaller than the electric dipole a\ for any finite 
value of e in the limit £ — > 0. For any non-zero value of 
the relaxation frequency 7 in ([5]) , the zero-frequency limit 
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FIG. 6: High-temperature Casimir free energy as a function 
of L/R calculated with the Drude model. The results do 
neither depend separately on L and R nor on the material 
parameters Xp and A 7 . 



FIG. 7: Thermal Casimir force correction computed with the 
Drude model. The conventions are the same as in Fig. [3] 
[Colors online]. 



of e is finite. In contrast with the perfect-reflector and 
plasma cases, the magnetic dipole contribution is then 
negligible 



LDrud 



lorn 4 

45c 



< a? rud 



(38) 



for any finite value of the Drude dc conductivity <tq = 
wp/j. The Drude sphere at low frequencies thus behaves 
as an inducible electric dipole for any finite value of <7o, 
corresponding to the Rayleigh scattering limit. 

The resulting high-temperature large-distance limit for 
the free energy reads 



T 



Drud 



hcR 3 

4\ T £ 3 



£> X T ,R. 



(39) 



Remarkably this result does not depend on the length 
scales Xp and A 7 characterizing the material response, 
whereas the corresponding plasma result (|3T))) clearly de- 
pends on Xp. One can show that this is always the case 
in the high-temperature limit At <S C, with J rDrud con- 
verging to the universal function of L/R shown in Fig. [51 
which is determined by the contribution of higher multi- 
poles i < £ max ~ R/L. 

The expression ((39)) corresponds to | of the value for 
perfect reflectors (f3"2")) , to be compared with the ratio \ 
found in the parallel-planes geometry [2J|, which results 
from the fact that the Fresnel coefficient tte vanishes at 
the zero frequency limit. Here the TE and TM contri- 
butions are redistributed into electric and magnetic mul- 
tipole spherical modes, thus explaining the change from 
\ to |. In fact, for perfect reflectors the magnetic dipole 
contribution proportional to b\ is one-third of the total 
free energy (|3"2"j). as can be surmised from (f2"T|) which 
shows that |6i| = |ai|/2. Since this contribution is negli- 
gible in the Drude model, the free energy is reduced by 
the factor |. 



Numerical results. An important consequence of the 
discussion presented above is that results from Drude 
and perfect reflector as well as plasma models are closer 
in the plane-sphere geometry than in the parallel-planes 
geometry. We have computed the Casimir force numer- 
ically for arbitrary (not too small) values of the surface 
distance L. For the plasma model, the thermal correction 
^pias (* n0 £ snown ) j s found to be close to the values for 
perfect reflectors shown in Fig. [3] On the other hand, the 
variation of $ Dmd calculated within the Drude model is 
remarkably different, as shown in Fig. [7J 

In contrast with the perfect reflectors and plasma 
model calculations, for which PFA always overestimates 
the thermal correction, PFA underestimates the thermal 
correction at short distances for the Drude model, while 
it overestimates it at long distances. The overestimation 
is, however, clearly smaller than for perfect reflectors. 
Since PFA results for plasma metals are above the val- 
ues calculated for Drude metals, the exact beyond-PFA 
Casimir force values at short distances for Drude and 
plasma models turn out to be much closer than predicted 
by the PFA-based theoretical models used in the analysis 
of experimental data. 

In order to highlight this striking feature, we plot in 
Fig. [5] the ratio between the thermal Casimir force val- 
ues F plas calculated with the plasma model and i^ Drud 
obtained with the Drude model, as a function of the dis- 
tance L, and for different values of the sphere radius R in- 
creasing from bottom to top. The fact that i^P las / j p Drud 
depends on R is a clear signature of the interplay be- 
tween geometry, temperature and dissipation. We also 
show the ratio computed within PFA (dashed line), which 
approaches 2 at large distances since the PFA result is 
proportional to the energy for parallel planes. On the 
other hand, the exact ratio approaches 3/2 for large radii 
R 3> Xp and more generally a value between 1 and 3/2 
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Appendix A: Matrix elements of the round-trip 
operator yVf (m) 



FIG. 8: Ratio between thermal Casimir force values calcu- 
lated with the plasma and the Drude models, as a function 
of surface separation L for different sphere radii. The solid 
curves from bottom to top correspond to increasing values of 
sphere radii. The PFA prediction corresponds to the dashed 
curve. [Colors online]. 



obtained from ([361 and (|39|): 

FP Us 3 / 1 cotha 



^Drud 



L> X T ,R. 



The factor 2 predicted by PFA in the large-distance limit 
is never approached by our exact results. 



VI. CONCLUSION 

In this paper, we have used the scattering approach to 
compute the Casimir free energy and force in the plane- 
sphere geometry, taking into account both the non zero 
temperature and the metallic nature of reflectors. 

For the simpler case of perfect metals at intermedi- 
ate distances, we observe from numerical computations a 
strong correlation between thermal and geometry effects, 
and negative entropy values for small spheres, which are 
clearly not related to dissipation but rather to geome- 
try itself. For small spheres, thermal photons provide 
a repulsive contribution, thus diminishing the total at- 
tractive Casimir force. Those results are endorsed by 
analytical derivations in the long distance regime. 

For the case of Drude metals, evidence of correlations 
between temperature and dissipation, qualitatively dif- 
ferent from those in the parallel plates geometry, is given. 
As a consequence, the results for the loss less plasma and 
full Drude models are closer to each other than in the 
parallel-plate geometry, with the long distance ratio of 2 
reduced to at most 3/2. If these results also hold in the 
experimental range R/L > 10 2 , it might diminish the dis- 
crepancy between experimental results and predictions of 
the thernral Casimir force using the Drude model. 



In this appendix, we derive explicit expressions for the 
matrix elements A^' m ^(^)i,2 given by (0. The different 
coefficients appearing in this equation are first calculated 
for real frequencies w. Given values of u> and k(fc, ip) de- 
fine two directions K(t9 , tp) in reciprocal space with 



ck 



sm( 



cost 



= ±- 



ck. 



k z = Woj 2 /c 2 — k 2 . 

( A1 ) 
We first consider the matrix elements implementing 

the change from the multipole to the plane wave basis. 

The free-space magnetic multipole fields have the Fourier 

representation [67] (K = KK.(0,ip))) 



r(M) 

Ldiril 



(K) 



,j(jf-^L^ 



L ic = l 



suit 



W/C y/Kjt + 1) 



\d v - (p k d e 



(A3) 



The electric multipole fields are in turn given by 

41(K) = -K xffi(K). 

In order to specify the plane-wave basis, we take the 
usual TE and TM polarizations adapted to treat the re- 
flection by the plane at z = 0: 



(?TE = z x k = tf> 
£tm = £te x K = 



(A4) 
(A5) 



The matrix element (k, — ,p\£, rn, P) are then obtained 
from the scalar product e p ■ ILaL(K ( _ )) with K'~) = 
(w/c)K(f9 _ , ip) after multiplication by the square root of 
the Jacobian dk z /d(uj /c) = u/(ck z ) : 



(k,-,TE\£mE) = 



(k,-,TE\£mM) = 



2inn^fuJ]c Y im {6 ,ip) 
kVK y/£(£ + 1) 
2ni dgYe m (6-,cp) 



(A6) 



(A7) 



^/k^Jfc ^/ £{£+!) 
(k, -, TM\£mE) = (k, -, TE\£mM) (A8) 

(k,-,TM\£mM) =-{*.,-, TE\imE). (A9) 

We now derive the Mie scattering matrix elements 
(£mP\7ts(0\^-! +iP)- We write the electric field in terms 
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of the Debye potentials (scalar fields satisfying Helmholtz 
equation) H( E ^(r) and I[( M )(r) for the electric and mag- 
netic multipole components. In order to have the De- 
bye potentials for a TE-polarized plane wave propagating 
along an arbitrary direction K(#, tp) (amplitude Eq), we 
take a rotation with Euler angles a = (p, ft = 6, 7 = 0. In 
terms of the coordinates x' , y' , z' corresponding to the ro- 
tated axis, the Debye potentials have the usual form [61] 
that corresponds to a plane wave propagating along the z' 
axis linearly polarized along the y' axis. To write in terms 
of the original coordinates, we use the matrix elements of 
finite rotations and find (ji(Kr) are the spherical Bcsscl 
functions 64]) 



En 



nI, TE (r,e,$) = ||£i 



Uir(2l+r 



Lm 



£(£ + !) 



■je(Kr) (A10) 



x e' im ^(dlA6) + d*(e))Y im (®,$) 



n 



K,T E (r,e,$) 



00 

On K 2-^1 



l,m 



>,/^W> (An, 



We use the same method to derive the Debye potentials 
for TM incident polarization: either we take the polar- 
ization along the x'-direction (instead of the y'-direction 
as done in the derivation for TE polarization) or else take 
the third Euler angle to be 7 = — ir/2 so that the rotated 
Oy' axis coincides with ctm instead of cte ■ This amounts 
to introduce an additional phase e Tn — ±i. Hence the 
potentials II K ' TM (r, 0,$) are obtained from (|A10j) and 

(TAUD by replacing < ±1 (0) -► (±»)<&,±i(0) 

Since the scattered field propagates outward from the 
sphere, the corresponding potentials are written in terms 
of the spherical Hankel functions h e (Kr) [64| . The De- 
bye potentials for the scattered field are then obtained 
by considering the boundary conditions at the surface of 
the sphere. In the resulting expression, each partial-wave 
term is multiplied by the corresponding Mie coefficient a,£ 
(electric multipoles) or bi (magnetic multipoles) |£2j. As 
expected from spherical symmetry, the Mie coefficients 
do neither depend on m nor on the direction of inci- 
dence. They are written in terms of the Riccatti-Besscl 
functions ipe((3) = f3ji((3), (i(f3) — f3h\ (8) evaluated at 
the size parameter (3 = ujR/c and at nf3 [61] : 



a f (uj) 
bt(ui) 






(A12) 



(A13) 



find the explicit multipole expansion for the scattered 
electric field (in position representation) (r|7^s|K, p), 
with p = TE, TM representing the incident polarization, 
which can also be cast into the formal decomposition 



(t\h s \k,p) = Y, 



ImP 



did 

2ttc 



r\ui£mP) (A14) 



x(uj£mP\n s \K, p), 

where the vector fields (r\u>£mP) are the inverse Fourier 
transforms of I w ^ m p(K) [see Eq. (|A2|) ]. Since the scatter- 
ing does not change the frequency, the matrix elements 
of IZs have the general form (K = k + k z z, k z > 0) 




[u£mP\n s \K, p) = 2ir8(K-u/c)d —(£mP\Tls\k,+, p) 

(A15) 
with the square root of the Jacobian dk z j d(yi j c) provid- 
ing once more the connection to our plane-wave basis 
|k, ±,p) u associated to a given frequency uj. 

By comparing the explicit expressions for (r|7?.,s|K, p) 
with the formal decomposition (|A14[1 and taking (IA15|) 
into account, we find 



(£mE\TZ s \K+,TE) 



lir(2i + l) 



ai{w) 



(A16) 



y k z uj/c 
x e-^(dlJ8+) + diJ9+)) 



(£mM\1ls\K+,TV) = iJ ^+^ bfa 



(A17) 



From the Debye potentials for the scattered field, we 



k z uj/c 
< e-^{dl A {9 + )-d^_ l {9 + )). 

For TM polarization, the matrix elements are obtained 
from (|A16j) and (|A17|) by substituting d l m±1 (9) -> 

(±»X ±1 W. 

The Casimir free energy may be written as an integral 
over the positive frequency semiaxis (which includes the 
evanescent sector < w < cfc). Using analyticity proper- 
ties of the reflection operators for plane and sphere over 
the upper complex frequency plane, we transform the in- 
tegral over real positive frequencies into an integral over 
complex frequencies u> = i£ + r], with £ running from +00 
to and 7/ — > + keeping the Matsubara poles i£ n outside 
the closed countour of integration employed in connec- 
tion with Cauchy theorem [20J . The resulting expression 
is given by ([6]), with the matrix elements M^ m \C)i.2 ob- 
tained by taking u> — > i£ and k z —> in, in (JA6I) - (1A9[) and 
(IA16[) - (|A17I) and plugging the results into ^\. The final 
explicit expressions are given in Sec. III. 
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